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ABSTRACT 



Context. MWC 297 is a nearby young massive B[e] star. The central star is attenuated by 8 magnitudes in the optical and has a high 

projected rotational velocity of 350 km s _I . Despite the wealth of published observations, the nature of this object and its circumstellar 

environment is not understood very well. 

Aims. With the present paper, we intend to shed light on the geometrical structure of the circumstellar matter that is responsible for 

the near- to mid-infrared flux excess. 

Methods. The H-band (1.6-2.0 /an), K-band (2.0-2.5 /an), and N-band (8-13 jim) brightness distribution of MWC 297 was probed 

with the ESO interferometric spectrographs AMBER and MIDI, mounted on the VLTI in Paranal, Chile. We obtained visibility 

measurements on 3 AMBER and 12 MIDI baselines, covering a wide range of spatial frequencies. Different models (parametrized 

circumstellar disks, a dusty halo) were invoked to fit the data, all of which fail to do so in a satisfying way. We approximated the 

brightness distribution in H, K, and N with a geometric model consisting of three Gaussian disks with different extents and brightness 

temperatures. This model can account for the entire near- to mid-IR emission of MWC 297. 

Results. The circumstellar matter around MWC 297 is resolved on all baselines. The near- and mid-IR emission, including the silicate 

emission at 10 micron, emanates from a very compact region (FWHM < 1.5 AU) around the central star. 

Conclusions. We argue that the extinction towards the MWC 297 star+disk system is interstellar and most likely due to remnants 

of the natal cloud from which MWC 297 was formed. Furthermore, we argue that the circumstellar matter in the MWC 297 system 

is organized in a circumstellar disk, seen under moderate (i < 40°) inclination. The disk displays no inner emission-free gap at 

the resolution of our interferometric observations. The low inclination of the disk implies that the already high projected rotational 

velocity of the star corresponds to an actual rotational velocity that exceeds the critical velocity of the star. This result shows that stars 

can obtain such high rotation rates at birth. We discuss the impact of this result in terms of the formation of high-mass stars and the 

main-sequence evolution of classical Be stars. 



Key words. Techniques: interferometric; circumstellar matter; Stars: pre-main-sequence; Stars: rotation; Stars: emission-line, Be; 
Stars: individual: MWC 297 

1. Introduction Il993t lHan et all fl998t ISkinner et al.l I1993I) . The image that 

emerges from these studies is far from uniform. 

At a distance of 250 pc and with a stellar mass of ap- _ , . . itt — — -r\ \ ,w& , i . 

, ,„ ,, } t^ . i l Nr.r.-r k » „, r/-i ™~, • r Based on an optical spectrum, Drew et al. ( 1997) have shown 

proximately 10 M dDrewetal. 1997h . MWC 297 is one of . . , ., ( \„,/„ m . , ■ , , 7T — — ", D1 -., , 

K , J . M r — • . ^ • . . convincingly that MWC 297 is a highly attenuated B 1.5 V star 

the closest massive young stars. It is therefore an ideal tar- ,. o\ tu a ■ a- » * \u e -><;n 

. , . , J . b . - __ (A v ~ 8). They derive a distance to the source of 250 pc. 

get for studying the formation process of massive stars. The i, , , , ... , , . ., 

° ■ , i , • T • • • Furthermore, these authors indicate that the star is a rapid ro- 

source has been the subject of many investigations, cover- »*/••• ocn cm -i\ a ^u » •. * *• 

, . , J » _. J b . tator (v sin i = 350 + 50 km s ) and suggest that its rotation axis 

ing the entire spectral range from X-rays to centimeter wave- .) , . . ., , c .i , ^ iroDTTMeru ai 

, . ^7^; — . ^^^r rV ~ : — '■ m ^»»^ t it; . ii^^-4 must be almost in the plane of the sky. The MERLIN 5 GHz data, 

lengths Vinketal. 2005; Damiam et al. 2006; Drew et al. 1997; , , . , F . , ,. / ., ,, , . 

a i n o T 11 ■< ^Ao » . ii i ™„„i i™ , ,i presented in the same article, display a north-south elongation, 

Andnllat & Jaschek 1998; Malbet et al. 2007; Teradaetal. F , DA ,. co D f XT p. ' / , 3 . . . ., . .,. ° 

along PA=165 E of N. Drew et al. point out that this elonga- 



2001; Benedettini et al. 2001; Acke & van den Ancker 200- 



Stl 



...... , , . , ,,,„, , . , -„,„- r , „ ,, ,, tion may reflect the equatorial plane of a disk-like circumstellar 

2006; Hennin g et all LL998; Manoi et al. 2007; iPalla & Prustil , , J , . ™/ on _ A .. .. . • . . .■ , , ., 

■ J structure around MWC 297. Both the high rotational velocity 

and the elongated radio structure are consistent with the pres- 
ence of an edge-on accretion disk, which is the byproduct of the 



Send offprint requests to: B. Acke stellar formation nf MWC 29? 

Based on observations made with ESO telescopes at the La Silla 
Paranal Observatory under program IDs 077.D-0071(B-C), 077.D- Near-IR interferometr ic observations, bot h broad-band 



0095(C-F), 079.C-0012(A-H) and 079.C-0207(A). (Mil lan-Gabet et al.l I200H lEisner et al. 2003; _|Monnier et al.l 



Postdoctoral Fellows of the Fund for Scientific Research, Flanders. 2006) and spectrally dispersed dMalbet et al.ll2007l) . have shown 
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that the radiation emitted by the hot circumstellar dust is con- 
fined to a very small region around the central star (ring diameter 
of 1.4 AU, M onnier et al . 2006). They suggest that the circum- 
stellar dust around MWC 297 is located in an accretion disk. In 
a non-accreting disk (i .e. a passive disk in the terminology of 
iDullemond et alJfeOOll) the temperature and density structure is 
predominantly determined by the stellar irradiation. Such a disk 
has a significant inner disk gap, due to dust sublimation close to 
the star. The source would display a much more extended, hence 
more resolved, brig htness distributio n on the sky. I n contr ast to 
the conclusion of iDrew et ail d 19971) . iMalbet et all J2007I) sug- 
gest an almost pole-on orientation of the accretion disk, based 
on the analysis of one spectro-interferometric AMBER measure- 
ment and regular spectroscopy of the Ha, H/3 and Bry lines. 
Their best model consists of a geometrically flat, optically thick 
accretion disk viewed under an inclination of approximately 20°. 
The hydrogen lines can be reproduced with a stellar-wind model 
which has a density enhancement close to the disk's surface. 
This could point to the presence of a disk wind. 

Spectropolarimetric observations of the Ha emission show 
no difference in linear polarization with respect to the continuum 
(i.e. line effect). This suggests a spherical or circle-symmetric 
emission region (lOudma iier & Drew 1999), which supports the 
low inclinati on found by Malb et et alj (120071) . 

Recently, Man oi et all d2007l) have presented 1.3 millimeter 
(mm) interferometric observations of MWC 297. They show that 
the mm spectral slope of the spectral energy distribution (SED) 
is indicative of grain growth in the circumstellar environment, or 
a very c ompact structure w hich is optically thick at mm wave- 
lengths. Manoi et al.l J20071) also study the CO emission around 
MWC 297. They find an extended optically thick CO cloud at 
vlsr ~ 10 km s which is resolved out in their interferometric 
observations. No compact CO emission region, associated to the 
mm-continuum of the MWC 297 system, is found. The authors 
interpret the latter as evidence of strong depletion of CO in the 
disk. 



2. Observations 

2.1. Interferometry 

MWC 297 was observed in the near- and mid-infrared with 
the inte rferometric instru ments AMBER dPetrov e t al. 2007J) and 
MIDI dLeinert et al.l I2003I) . Both instruments are mounted on 
the VLTy. The observations were performed with the 1.8m 
Auxiliary Telescopes (ATs) in different baseline configurations. 
In this paper, we use the VLTI nomenclature to refer to the dif- 
ferent telescope stations and baseline settings. 

AMBER simultaneously combines the light of three tele- 
scopes and currently operates in H- and K-band (1.6-2.5 /mi). 
The fringes are spectrally dispersed. Our data were taken in the 
low spectral resolution mode (R~30, LR-HK mode). The data 
re duction was perform ed according to the methods described 
in iTatulli etalj d2007l) . Each observation consists of a number 
of frames, i.e. exposures with short integration times. We have 
tested different frame-selection algorithms (e.g. the suppression 
of frames with low photon counts in the photometric channels; 
the selection of frames with statistically significant piston es- 
timates; the selection of frames with high signal-to-noise visi- 
bility estimates), which were applied to both the target and the 
calibrator star. Typically, 60-100% of the frames were selected. 
Little difference is seen between the resulting calibrated vis- 
ibilities, which confirms the excellent quality of the AMBER 

1 http://www.eso.org/projects/vlti/ 



data. Furthermore, we have repeated the AMBER observations 
on the same baselines to check their consistency. The results of 
both -independent- observations are equal within the error bars. 
We have averaged both measurements to obtain the final set of 
AMBER visibility data. 

Thanks to the combination of 3 telescopes, AMBER also 
provides a spectrally dispersed closure phase (CP) measure- 
ment. To estimate the uncertainty on the latter, we take the stan- 
dard deviation of the closure phase measurement of the calibra- 
tor, which is expected to have a zero closure phase. The clo- 
sure phase of MWC 297 is within the error bars equal to zero 
(|CP| < 5°). Brightness distributions on the sky that display devi- 
ations from centro-symmetry should always display a non-zero 
closure phase, except at specific combinations of positions in 
the (u,v)-plane where phases happen to sum up to zero. This is 
also true if the baselines in the triangle are aligned, which is 
the case for our observations. Hence, the zero AMBER closure 
phase strongly suggests the circumstellar geometry of MWC 297 
is point-symmetric along the position angle of the baselines, at 
the spatial resolution of our observations. 

We have extracted the H- and K-band AMBER spectrum 
of MWC 297 from the three photometric channels. The flux 
in each channel is summed over all frames, and calibrated us- 
ing the total flux in the corresponding channel of the calibra- 
tor star. A photospheric model of the calibrator star is applied 
to convert to flux units. Finally, the absolute flux level of the 
AMBER spectrum is determined by multiplicative scaling to the 
mean value of the Infrared Space Obser vatory Short Wavelength 
Spectrograph (ISO-SWS) s pectrum dBenedettini et alj 1200 it 
lAcke & van den Anckerl 2004) in the overlap region between 
both (2.36-2.5 /mi). The error on the spectral shape of the spec- 
trum is determined from the rms of the three photometric chan- 
nels. We note that the aperture of the ISO-SWS instrument at the 
shortest wavelengths is 14" x 20", much larger than the seeing- 
limited resolution of a 1.8m telescope. However, ground-based 
photometric measurements are in agreement with the ISO flux 
levels. This indicates that the fraction of extended emission on 
a large scale is negligible, and that the scaling of the AMBER 
spectrum to the ISO spectrum is allowed. 

MIDI is a two-beam combiner which operates in N-band (8- 
13 /mi). The fringes are spectrally dispersed at a resolution of 
30 (HIGH-SENS PRISM mode). We have used the MIA+EWS3 
data reduction packages to obtain calibrated visibilties. The 
so-called incoherent (based on power-spectrum analysis of the 
fringes) and coherent (based on the information contained in the 
spectrally dispersed fringes) data reduction results are consis- 
tent with each other, which indicates the good quality of the 
data set. Furthermore, the coherent analysis provides differen- 
tial phases. The latter are within the error equal to zero on all 
MIDI baselines. This suggests that also the N-band emission re- 
gion, including the silicate feature, does not deviate much from 
centro-symmetry. 

The MIDI spectrum is determined in the following manner: 
The raw spectrum of the science target is obtained by adding all 
photometry files. It is flux-calibrated as described bv lvan Boekell 
(2004) with the corresponding quantity of two calibrators, or two 
separate measurements of the same calibrator, observed in the 
same night as MWC 297 and spaced in airmass. Again, model 
photospheres were applied to compute the theoretical spectrum 
of the calibrators. The angular extent of the calibrators was used 
to scale the spectra to the right absolute flux level. This proce- 
dure was obviously only possible in nights were more than one 



http://www.strw.leidenuniv.nl/~koehler/MIA+EWS-Manual/ 
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Fig. 1. Interferometric data. Left, top to bottom: the calibrated AMBER visibilities, closure phases and spectrum. Right, top to 
bottom: the calibrated MIDI visibilities, differential phases and spectrum. Different colors refer to different baselines. In the top left 
panel, the AMBER data corrected for the flux contributio n of the central star are overp lotted in grey (see Sect. 13. 11 1. The grey line 
in the bottom panels is the ISO spectrum of MWC 297 (lAcke & van den Ancker| [2004). Representative error bars are overplotted, 
except for the MIDI differential phases; the average uncertainty on this quantity is 1 degree, consistent with the spread of the data 
around zero. [This figure is available in color in the electronic version.] 



calibrator was observed in the same MIDI setting. We obtain 
9 independently calibrated spectra. The final MIDI spectrum is 
the average of these, with the error being the standard deviation. 
It is striking that the independently calibrated ISO spectrum and 
MIDI spectrum of MWC 297 have the same shape and flux level, 
consistent with the IRAS 12 pm flux, despite the very different 
aperture size or field-of-view of the respective instruments. This 
argues in favor of a mid-IR emission region which is spatially 



unresolved by a single telescope. The entire N-band emission 
region is hence sampled by the MIDI interferometer. This con- 
clusion is confirmed by a 11.85 pm VISIR image of MWC 297, 
in which the source is resolved, but on a much smaller scale than 
observable with a 1.8m telescope (see Sect. 14.2b . 

The log of the interferometric observations can be found in 
Table [TJ The calibrator stars are indicated as well. The result- 
ing visibilities, closure phases and spectra are plotted in Fig. Q] 
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The visibilities are shown as a function of spatial frequency B/A, 
with B the projected baseline length and A the wavelength of the 
observations. The spatial frequency is a measure for the resolv- 
ing power obtained with a certain measurement. To guide the 
reader who is not an interferometry expert: each measurement 
at a certain baseline consists of a spectrum of visibilities and 
phases, which not only sample a different wavelength, but also 
a different spatial frequency and direction on the sky. Visibility 
measurements at the same spatial frequency may therefore dif- 
fer due to a different baseline position angle and/or the spectral 
dependence of the source's emission on the sky. In the top right 
panel of Fig. Q] the MIDI visibilities are presented. Here, it can 
be seen that the 10 micron silicate feature leaves a spectral im- 
print on the visibilities. We will come back to this in Sect. 13.2.31 



2.2. Additional data 

We have expanded our data set with previously published pho- 
tometric and spectroscopic data. In Fig. [7] the dereddened 
SED is shown, includ ing the ISO spectrum. The photometri c 
data were taken from iBerrilli et ail dl987l) . iManningsl dl994l) . 
Ide Winter etail d2001l) , and the 2MASS and IRAS catalogues. 
The central star is modeled with a Kurucz model for a B1.5V 
star (T e ff = 25,400K, logg = 4.0). The stellar radius, obtained 
from the fit of the photosphere model to the UV and optical 
photometry, is 8 ±2 R . The R- and I-band photometric points 
were excluded from this fit, since excess flux is present at these 
wavelengths (~40% and 20% respectively). The excess in R can 
be entirely attributed to the extremely strong He emission in 
MWC 297. Also the I-band excess is due to hydrogen emis- 
sion li nes. The Paschen lines measured by I Andrillat & Jaschekl 
( 1998) account for the ent ire excess. We hav e applied the inter- 
stellar dereddening law of lSavage & Ma this ( 1979), e xtended to 
infrared wav elengths with the law described by Steenman & The 
(fl989l[T99Th . We set R y = 3.1 andE(B-V ) = 2.5, leading to A v 
= 7.75, consistent wifh lDrewetaH d 19971) . The AMBER, MIDI 
and ISO spectra are plotted in Fig. [2] After dereddening, the ISO 
spectrum displays a weak but significant 10 micron silicate emis- 
sion feature. The infrared excess luminosity Lir is only 5% of 
the stellar luminosity L*. 

In Section l4~Tl we argue that the extinction of the central star 
is most likely due to a foreground interstellar molecular cloud, in 
between the observer and the MWC 297 system. Therefore, the 
real spectrum of the MWC 297 system is the dereddened spec- 
trum. In the following, we use the dereddened AMBER, MIDI 
and ISO spectra, unless otherwise specified. We stress that the 
dereddening has no direct consequence for the interferometric 
measurements, since the latter are only sensitive to the normal- 
ized brightness distribution of the source at each wavelength. 

Diffraction-limited images of MWC 297 in a broad-band fil- 
ter centered at 11.9 yum, and intermediate- width filters centered 
on 11.3 and 11.9 /mi were obtained with t he VLT Imager and 
Spectrometer for the mid-Infrared (VISIR, lLagage et al] r2004) 
during the night of April 13-14, 2007. Subtraction of the ther- 
mal emission from the sky, as well as the telescope itself, was 
achieved by chop-nodding in the North-South direction with a 
throw of 30". The effective integration time was 90 seconds in 
each filter. Calibrator measurements in the same settings were 
obtained shortly before the science observation. All observa- 
tions were obtained when the sources were close to zenit (air- 
mass ~ 1.05). The full-width-at-half-maximum(FWHM) of the 
calibrators' broad-band image is 0.33±0.01", the FWHM of the 
MWC 297 image is 0.41". 
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Fig. 2. ISO (black line), AMBER (dots) and MIDI (dots) spec- 
tra of MWC 297. The AMBER spectrum was scaled to the ISO 
spectrum in the overlap region between both spectra. The MIDI 
spectrum has been calibrated independently and agrees well with 
the ISO spectrum. In grey, the ISO spectrum is shown, dered- 
dened with the same extinction law as the optical photometry. 
A mild but significant silicate emission feature is present at 10 
micron. 



MWC 297 has been observed during a Science Verification 
run of the Atacama Pathfinder Experiment (APEX0). The 
APEX-2A heterodyne receiveiQ mounted on this telescope pro- 
vided a raster map of the CO (3 - 2) emission line at 345.79 
GHz in a 1' x 1 'field around the source, with a beam size of 
17.5". No spatial dependence of the spectrum has been detected 
in the field, indicating that the detected emission lines are due to 
a CO structure, significantly large r than the mapped area. This 
is in agreement with the results of Ma noi et all d2007l) . who at- 
tribute this emission to interstellar clouds. In Fig. [3] the field- 
averaged spectrum is shown. No features are detected at high 
red- or blueshift. The latter excludes the presence of a strong 
CO outflow i n MW C 297, on top of the non-detection of CO by 
Manoi etall d2007l) in the circumstellar disk. 



3. Analysis 

3.1. Simple geometric modeling 

We have fitted simple analytic models to the visibility data. 
Because the AMBER closure phases and the MIDI differen- 
tial phases are consistent with zero on all baseline settings, we 
have only applied point-symmetric and concentric models. The 
brightness distributions in the H-, K-, and N-bands may be domi- 
nated by different components. Therefore, we have first analyzed 
the data in the near- and mid-IR separately. In the last part of this 
section, we confront the results and construct a model which fits 
the interferometric observations and spectra in both wavelength 
ranges simultaneously. 



3 This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between the 
Max-Planck-Institut fur Radioastronomie, the ESO, and the Onsala 
Space Observatory. 

4 http://www.apex-telescope.org/heterodyne/het345/ 
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Table 1. The log of the interferometric observations, sorted according to interferometric instrument and increasing baseline length. 
The first seven columns refer to the observations of MWC 297, the last four to the calibrator's properties. As a reference, the 
wavelength-averaged calibrated visibility <V>a is given. The flux of the calibrator is the 2MASS Ks-band magnitude for AMBER 
observations and the IRAS 12 /jm flux for MIDI. PB = projected baseline length; PA = projected baseline angle. 



Night 


UT 


Instrument 


Setting 


PB 


PA 


Airmass 


<v>,. 


Calibrator 


Sp.T. 


Diameter 


Flux 


yyyy-mm-dd 


hh:mm:ss 






m 


° EofN 










mas 


mag or Jy 


2007-04-14 


08:09:13 


AMBER 


E0-G0 


14.2 


68.4 


1.16 


0.78 


HD 166460 


K2III 


1.40 


2.58 


2007-06-19 


07:10:31 


AMBER 


E0-G0 


14.9 


73.0 


1.20 


0.86 


HD 166460 


K2III 


1.40 


2.58 


2007-04-14 


08:09:13 


AMBER 


G0-H0 


28.5 


68.4 


1.16 


0.60 


HD 166460 


K2III 


1.40 


2.58 


2007-06-19 


07:10:31 


AMBER 


G0-H0 


29.9 


73.0 


1.20 


0.59 


HD 166460 


K2III 


1.40 


2.58 


2007-04-14 


08:09:13 


AMBER 


E0-H0 


42.7 


68.4 


1.16 


0.46 


HD 166460 


K2III 


1.40 


2.58 


2007-06-19 


07:10:31 


AMBER 


E0-H0 


44.4 


73.0 


1.20 


0.41 


HD 166460 


K2III 


1.40 


2.58 


2007-04-12 


08:22:50 


MIDI 


E0-G0 


14.3 


68.5 


1.15 


0.79 


HD 148478 


M1.5I 


25.92 


3200 


2006-04-20 


06:02:07 


MIDI 


D0-G0 


20.1 


56.3 


1.58 


0.63 


HD 146051 


M0.5III 


9.83 


150 


2007-04-17 


06:50:30 


MIDI 


G0-H0 


23.4 


61.9 


1.37 


0.60^, 


HD 82668 
HD 187642 


K5III 

A7V 


6.99 

3.22 


73 
33 


2006-04-20 


07:57:03 


MIDI 


D0-G0 


28.9 


68.8 


1.13 


0.51 


HD 168454 


K3III 


5.78 


62 


2006-04-21 


09:09:21 


MIDI 


D0-G0 


31.8 


72.1 


1.07 


0.46 


HD 167618 


M3.5III 


11.33 


214 


2006-04-22 


09:59:19 


MIDI 


D0-G0 


31.9 


73.1 


1.09 


0.38 


HD 168454 


K3III 


5.78 


62 


2007-04-19 


08:00:32 


MIDI 


E0-H0 


43.3 


68.8 


1.14 


0.30 


HD 187642 


A7V 


3.22 


33 


2006-05-23 


05:06:38 


MIDI 


A0-G0 


52.6 


66.0 


1.23 


0.27 


HD 167618 


M3.5III 


11.33 


214 


2006-05-25 


05:12:08 


MIDI 


A0-G0 


54.5 


67.1 


1.19 


0.25 


HD 146051 


M0.5III 


9.83 


150 


2006-05-25 
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MIDI 


A0-G0 


64.0 


72.7 


1.07 


0.18 


HD 146051 


M0.5III 


9.83 


150 


2007-05-10 


05:45:56 


MIDI 


G1-H0 


66.8 


172.9 


1.27 


0.14 


HD 167618 


M3.5III 


11.33 


214 


2007-05-08 


06:41:00 


MIDI 


D0-G1 


71.5 


129.5 


1.14 


0.13 


HD 167618 


M3.5III 


11.33 


214 



" The calibrators HD 82668 and HD 187642 were observed at high airmass (~2), which hampers a reliable calibration. We have averaged the 
calibrated visibilities. The large errors (Fig.fTJ were deduced from the standard deviation. 
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Fig. 3. APEX spectrum of the CO (3 - 2) emission line, aver- 
aged over a 1' x 1' field around MWC 297. The three detected 
features in the VLSR-range between 2 and 13 km s -1 are spatially 
extended and most likely produced in interstellar clouds. No fea- 
tures exceeding the noise level are detected at any other velocity. 



3.1.1. Fit to the MIDI data 

Most of the MIDI data are obtained on baselines along the VLTI 
0-line (i.e. projected baseline angle ~70°). We are hence pre- 
dominantly sensitive to the N-band brightness distribution pro- 
jected onto this direction. The flux contribution of the central star 
in the N-band is expected to be less than 0.5% and is therefore 
neglected in the fits described below. 

For each observational wavelength, the MIDI visibility curve 
shows a steady decline with increasing spatial frequency. Strong 
lobes are absent, although a drop-off may be present at spatial 



frequencies around 30 cycles per arcsec. One-component uni- 
form disk (UD), uniform ring (UR) and Gaussian disk (GD) 
models cannot fit this shape. Two-concentric -component mod- 
els (UD + point source, GD + point sourcq^, two UDs, two 
GDs) do better. In general, the GD models do better than the UD 
models (a factor two in reduced chi-square). In the best-fit ring 
models, the inner gap is always smaller than 1 mas (0.25 AU 
at 250 pc). The ring models are thus essentially equal to the 
uniform-disk models. We have summarized the best-fit param- 
eters of some illustrative models in Table [2] The amount of free 
parameters in the two-component models was alternated. First, 
we have fitted the extent of both components, and allowed for a 
wavelength-dependent linear slope of the flux ratio (model type 
A in Tabled. The extent of the components was constrained to 
be the same over the N-band. Hence, this model is fully deter- 
mined with only 4 parameters. Second, the extent of the com- 
ponents is again wavelength-independent, but now the flux ratio 
at each wavelength is a free fitting parameter (model type B). 
The total number of parameters is 2 plus the number of wave- 
length points. Finally, we have also fitted models in which both 
the extent and the flux ratio are wavelength-dependent (model 
type C). Models including point sources have worse x^ ed values 
and are therefore excluded from the table. The best model is a 
2GD model with all parameters /1-dependent, with a^Lj close to 
1. The best 2UD model overfits the data (yjL = 0.6). 

A number of authors (e.g. iMalbet et alTl2007b iManoi et all 
2007) have suggested that MWC 297 is surrounded by an accre- 
tion disk. Therefore, we have fitted an alpha-disk (AD) model, in 
which it is assumed that each annulus around the central star ra- 
diates as a blackbody. The temperature at each distance from the 
star is parametrized with a power law T(R)ocR°\ The six free pa- 
rameters of this model are the inner and outer radius, the temper- 



5 Note that the point source in these models does not represent the 
central star, but an unresolved disk component. 
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Table 2. Results of some illustrative analytic model fits to the 
MIDI visibilities. The extent of the model components CI and 
C2 is the wavelength-averaged FWHM for GDs and diameter for 
UDs. The wavelength-averaged flux ratio is included as well. #p 
is the number of fitting parameters in the model, )Q, the reduced 
chi-square value. #A is the total number of wavelength bins in 
one MIDI observation. See text for an explanation on the model 
types. 



Model 


Best-fit parameters 


#p 




type 


Extent C 1 Extent C2 Flux ratio 
mas mas C2/C1 







2 UDs A 


65 


2 GDs A 


38 


2UDsB 


65 


2GDsB 


40 


2UDsC 


96 


2GDsC 


79 



23 
10 
23 
12 
32 
20 



0.74 


4 


4.2 


0.28 


4 


2.8 


0.79 


#A+2 


3.6 


0.38 


#/t+2 


1.9 


1.44 


3x« 


0.6 


1.35 


3x#A 


1.1 



Table 3. Results of alpha-disk (AD) fits to the MIDI visibilities. 
The alpha disk is determined by the inner and outer disk radius 
R; n and R out , the disk inclination i and position angle PA, and 
the temperature law T = To (R/Ri n ) ff . #p is the number of fitted 
parameters in the model, x 2 kA the reduced chi-square value of the 



visibilities. 



Model 




Best-fit parameters 


#p 


2 


AD 


To = 


1510K 


a = -0.68 


5 


4.9 




Rin ' 


= 4 mas 


i = 40° 








AV>llt 


= 320 mas" 


PA = 160° 






Standard 


To = 


1650 K 


a = -0.75" 


4 


6.9 


accretion 


Rin = 


= 4 mas 


i = 45° 






disk 


K-out 


= 320 mas" 


PA = 150° 







Table 4. Results of some illustrative analytic model fits to the 
AMBER visibilities, corrected for the flux contribution of the 
central star. The nomenclature is the same as in Table|2] See text 
for details. 

Model Best-fit parameters #p ~X^ A 

type Extent C 1 Extent C2 Flux ratio 

mas mas C2/C 1 



UD+P 


14.6 


UR+P 


0.2 


GD+P 


8.3 


2 UDs 


23.2 


2 GDs 


14 



0.77 


2 


2.5 


0.44 


3 


2.5 


0.46 


2 


1.6 


3.6 


3 


0.9 


2.8 


3 


0.7 



7.3 

7.3 
4.3 



3.1 .2. Fit to the AMBER data 

Our observations in the H- and K-band are sensitive to the hottest 
circumstellar matter present in the MWC 297 system. At these 
wavelengths, also the central star is a contributor to the total flux. 
The stellar radius of 8 R corresponds to an angular radius of 
0. 15 mas. At the spatial resolution of our AMBER observations, 
the star is therefore completely unresolved and can be treated as 
a point source. 

Supported by the zero closure phase, we assume that the pho- 
tocenter of the circumstellar emission coincides with the central 
star. For each wavelength bin, we can then extract the visibil- 
ity V csm of the circumstellar matter from the measured visibility 

Vmeas Using 



Fixed a priori. 



v„ 



F+Vt + F^Vt 



csm y csm 



F+ + F c: 



ature at the inner radius and the power of the temperature law, the 
inclination of the disk and the position angle of the major disk 
axis on the sky. In a standard accretion disk, the temperature 
drops off with radius proportional to R" 75 . An accretion disk 
can therefore be modeled with an alpha disk in which the power 
is set to -0.75. The parameters of the best-fit AD models to the 
MIDI data are summarized in Table [3] The outer radius of the 
disk is ill-defined by our data set, i.e. equally good fits can be ob- 
tained with any value for this parameter. We have th erefore fixed 
R out to the maximum outer radius determined by iManoi et alj 
(2007, 80 AU) from mm-observations. The AD models produce 
a spectrum, which is constrained by the observed MIDI spec- 
trum. Reasonable fits can be obtained, although the simple two- 
component models presented in Table [2] appear to be much bet- 
ter. 

The main conclusion in this section is that simple analytic 
two-component models do a much better job in reproducing 
the observed MIDI visibilities than single-component models, 
a uniform ring or an alpha disk. Two circumstellar structures 
are present in the MWC 297 system; a smaller structure with 
a radius of ~10 mas (2.5 AU at 250 pc) and a larger structure 
of ~40 mas (10 AU). Note that both components are extremely 
compact, given the strong radiation of the central B 1 .5 V star. In 
fact, both structures reside entirely within the 1500K dust sub- 
limation radius of 15 AU for small (0.1 yum) amorphous olivine 
grains. 



with V+ = 1 , and F* and F csm the (uncorrected) flux of the star 
and the circumstellar matter respectively. The latter are com- 
puted from the AMBER spectrum and the photospheric model 
of the central star that fits the SED. In H-band, the stellar flux 
contribution adds up to 20%, in K to 10%. Hereafter, we refer to 
y csm as the AMBER visibilities, unless otherwise indicated. 

Simple analytic models were fitted to the AMBER data. In 
Table|4]we summarize some illustrative results. All models have 
wavelength-independent parameters, apart from the initial cor- 
rection of the visibilities for the star-to-disk flux ratio. Increasing 
the complexity of the models, e.g. by allowing a flux ratio of 
the components that varies with wavelength, only decreases the 
quality of the fit in terms of reduced chi-square due to the small 
number of AMBER data points. We note that the best-fit UR 
models prefer a disk which has no inner gap (Rj n < 0.2 mas » 
R*). As a result of this, the UR and UD fits are basically the 
same, with comparable %?, values. It is remarkable that two- 
component models, with a fixed extent and flux ratio over the H- 
and K-band, explain the AMBER data very well within the ob- 
servational errors. A tit with an alpha-disk model to the AMBER 
data alone is not possible, given the small number of AMBER 
points, and the fact that all AMBER baselines are aligned: the 
inclination and position angle of the disk are not constrained, 
the other parameters vary over a broad range in parameter space. 
The parameter a, however, can be restricted based on the slope 
of the visibility curve with wavelength. It appears to lie prefer- 
entially in the range between -0.6 and -0.5, instead of around 
the standard value for an accretion disk (-0.75). 
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3.2. Combining AMBER and MIDI 

3.2.1. Disk models 

The major advantage of the combination of near-IR AMBER 
and mid-IR MIDI data is the sensitivity to the temperature gra- 
dient of the circumstellar matter. Analytic alpha-disk models can 
therefore be subjected to a constraining test. The models are not 
only constrained by the visibilities, but also by the spectra in H, 
K and N. The best-fit standard accretion disk (a fixed to -0.75) 
has a reduced chi-square of more than 60 and can therefore be 
ruled out. Leaving the power of the temperature law free, the best 
fit is obtained with R in = 0.42+0.06 AU = 1.7+0.2 mas, R out = 
20+20 AU = 80±80 mas, a temperature at R in of 1600+250K 
and a power a=-0. 55+0.03. The inclination and position angle 
of the best-fit alpha disk are (=20+20° and PA= 180+50° respec- 
tively. The latter, as well as the outer radius, are ill-defined by the 
data set. The reduced chi-square of this model is 8.8. The uncer- 
tainties on the best fit parameters are derived from fits to virtual 
data sets, Monte-Carlo simulated around the original data. 

Alpha-disk models fail to reproduce the observations, be- 
cause the presence of an inner gap is needed in the models to 
keep the temperatures in the inner disk physical (i.e. not exceed- 
ing a few 10,000K). The visibilities show however that no gap 
in the brightness distribution is present at the resolution of our 
observations. No secondary lobes are observed in the visibility 
curves, which is e.g. witnessed by the tendency of the uniform 
ring model fits to have an inner radius close to the stellar surface. 
Also passive disk models with an inner disk gap (Dullemond 
2002, and subsequent papers), in which the central star is the 
only heating source of the circumstellar matter, can a priori not 
explain the observations. If a hypothetical passive disk around 
MWC 297 is seen under a moderate inclination (0° < i < 70°), 
the inner rim at dust sublimation temperature would dominate 
the near- to mid-IR brightness distribution. Given the large sub- 
limation radius, a detectable gap would show up. If the system 
is seen edge-on, its appearance on the sky changes. Although 
present, a physical gap in the disk would not be detected in the 
visibilities. 



3.2.2. A dusty outflow or halo? 

It is suggested that the spectral characteristics of young stellar 
objects can sometimes be explained without the need for a 
disk-like geome try but instead a ssuming a spherical distribution 
of the dust (e.g. lVinkovic et alj|200 6). Given the unclear origin 
of the circumstellar material around MWC 297, we deem it nec- 
essary to test this hypot hesis. To this end, we used the radiative 
transfer code MoDust (iBouwman et al.ll2000t lBouwmanll2001l) 
with optical constants from the Jena database^. 

A dusty outflow 

Mass loss through a dusty wind such as those found around 
evolved stars is an a priori improbable source of the IR excess, 
as the observed features are crystalline rather than amorphous, 
with a large continuum contribution suggesting the presence of 
at least micron-sized grains. Our attempts to model SED and 
visibilities simultaneously with an outflow model failed in all 
possible approaches, regardless of the dust species and grain 
characteristics used: (1) An outflow model matching the SED 



characteristic^] can be found, with an inner dust radius compat- 
ible with its condensation temperature (R; n = 907?* = 3.5 AU, 
Ti n = 1500 K), but the derived visibilities are orders of mag- 
nitude too small at all wavelengths. (2) Focussing instead on 
the observed spatial dimensions of the circumstellar structure, 
it is possible to reproduce the visibilities in the mid-IR in 2 
ways: either with a small inner dust radius or starting from 
the SED-matching model and reducing the outflow to a thin 
shell at the inner radius. In the first case, the dust which is 
unshielded from the stellar radiation reaches unrealistic temper- 
atures of 3000 K and, moreover, this hot dust generates such a 
strong near-IR excess that the AMBER visibilities can not be 
reproduced. This is the case even when using only mm-sized 
or oxide grains which resist heating much more than smaller 
and/or silicate grains. In the second case, the SED is no longer 
matched, and the dust which is at condensation temperature 
is still too bright and/or resolved in the near-IR. (3) No way 
was found to reproduce the AMBER visibilities, even without 
the constraint of the SED and MIDI visibilities, with anything 
remotely resembling a dusty outflow. 

Other density distributions 

We attempted fits with other density distributions as well, 
using power laws (p(R) = po R? with p the dust density, R the 
distance from the central star and /3 the free parameter determin- 
ing the density gradient) with/? ranging from -3 to +3. Again, the 
inner and outer radius were varied but all to no avail: requiring 
radiative equilibrium for the dust makes it impossible to have it 
at a temperature of only 500 K (the brightness temperature of the 
larger component, and the temperature preferred for that compo- 
nent by the analytical models presented in Sect. 13.2.31 ) so close 
to the central star without some form of effective shielding. All 
fitted models have reduced chi-square values over 100. 

We conclude from these attempts that a dusty halo, be it an 
outflow or some other density distribution, can explain none of 
the components observed with AMBER and MIDI, in spite of 
good spectral agreement. 

3.2.3. Approximating the brightness distribution 

No tested physical model (accretion disk, halo, passive disk) is 
able to explain the observations. Therefore, we aim at approxi- 
mating the near- to mid-IR brightness distribution of the source 
using a combination of analytic models. This approach focusses 
on the fitting of all interferometric data as a whole and, as a by- 
product, delivers an image of the source close to its real appear- 
ance on the sky. One should however keep in mind that the ana- 
lytic components used in the fits do not necessarily have physi- 
cal meaning. They are merely building blocks for our image. We 
will come back to this in Sect. [4] 

Two component models fit the AMBER and MIDI data well 
within the errors. The AMBER visibilities, corrected for the con- 
tribution of the unresolved central star of the system, are flux- 
dominated by a small component. In the MIDI data, on the other 
hand, the most extended component is dominant. The extent 
of the largest component in the AMBER fits, and the smallest 
component in the MIDI fits, appears to be of the same order. 
Furthermore, the closure and differential phases, in HK- and N- 
band respectively, suggest a point-symmetric source geometry. 
Inspired by the AMBER and MIDI fits, we have built a concen- 



http://www.astro.uni-jena.de/Laboratory/Database/databases.html 



7 Due to its large inner radius, this model hardly obscures the photo- 
sphere, showing that a low circumstellar Ay alone does not necessarly 
require a disk-like dust distribution. 
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Table 5. The parameters of the best-fit three-GD model. Top: Fit 
to all interferometric data simultaneously. Middle: Fit to subset 
only. Bottom: Inclination and position angle (PA) of the disk's 
major axis, based on the visibility measurements obtained on 
the two baselines with different baseline position angles. The 
reduced chi-square ^ 2 ed for each fit is given. See text for details. 
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trie three-component model, consisting of three Gaussian disks 
(3GD), to fit all visibilities simultaneously. We assume that two 
of these components are hot enough to radiate in the near-IR, 
and that the third one represents the cooler material, which is 
responsible for the 10 micron silicate feature and does not con- 
tribute to the AMBER flux. Since the brightness temperatures 
of the two HK-band components is high (as expected from the 
AMBER fits), the N-band flux contribution of the two hot com- 
ponents in the 3GD model is represented by the Rayleigh- Jeans 
tail of a blackbody (i.e. F„ oc A' 2 ). This approximation is made 
to suppress the number of model parameters. In summary, the 
model has 7 free parameters. The FWHMs of the GDs consti- 
tute three of those. Over the AMBER wavelength range (1.6— 
2.5 /im), we allow for a linear slope of the flux ratio of the two 
hottest components, implying two additional model parameters. 
The last two parameters control the absolute flux levels of the 
blackbody tails of the two hottest components at 10 micron. The 
flux contribution of the coolest component is then fixed to the 
N-band spectrum minus these blackbody contributions. At this 
stage, we do not allow for asymmetry due to inclination. This 
will be discussed later in this Section. 

We fit this model to all visibility data simultaneously. The 
best-fit parameters and errors are summarized in the upper part 
of Table [5] The three GD components were labeled A, B and C 
from small (hot) to large (cold). The parameters are the FWHM 
of each component, the linear flux ratio (FR) of component A 
and B in the near-IR, and the absolute flux levels (FL) of compo- 
nents A and B at 10 fim. The uncertainties on the best fit parame- 
ters are derived from fits to one hundred Monte-Carlo simulated 
data sets, around the original data. 

In our data set, two baselines have significantly different po- 
sition angles than the rest (PA=170° and 130° versus PA«70° 
respectively). If present, strong asymmetries in the brightness 
distribution should be observable at these baselines. Therefore, 
we have repeated the fit on a smaller data set (hereafter called 



the subset of data), limited to the data obtained at PA»:70°. 
Also these best-fit parameters are included in the middle part 
of Table[5] We note that the reduced chi-square is noticably bet- 
ter (y 2 ed = 4.1 versus 4.5, respectively), although only 85% of 
the data was used, and the number of free model parameters re- 
mained the same. The extent of the three Gaussian disk compo- 
nents are comparable in the two fits, and are 4, 12 and 41 mas 
respectively. This corresponds to 1, 3 and 10 AU at a distance of 
250 pc. 

The blackbody brightness (BB) temperatures of the compo- 
nents are deduced from their relative flux contributions to the 
AMBER+MIDI spectrum. The errors are determined from the 
rather large errors on the absolute flux level on the spectra, and 
the uncertainties on the model parameters. The smallest com- 
ponent has a BB temperature of 1700±200K, the middle com- 
ponent 930±80K. As a consequence of our modeling approach, 
the spectrum of the largest component contains the 10 micron 
silicate emission feature, as well as part of the underlying con- 
tinuum. We have determined the temperature of this component, 
assuming that its flux is exclusively produced by 0. 1 //m olivine 
grains. The best fit was obtained with a temperature of 520+20K 
and a column density larger than 0.02 gcirT 2 . The latter value, 
in combination with the spatial extent of the component, corre- 
sponds to a mass in small olivine grains of at least 2 x 10~ 7 M Q . 

To get a handle on the possible elongation of the source's 
brightness distribution - e.g. due to inclination - we have taken 
the 3GD model that fits the subset of data, and applied it to the 
visibility data obtained on the two baselines with different po- 
sition angles. We insert one additional fitting parameter, a mul- 
tiplicative factor that can change the extent of the three compo- 
nents simultaneously depending on baseline position angle. All 
relative quantities (FWHM ratios, flux ratios) were fixed to the 
best-fit parameters, while the absolute scaling of the extent of 
the model was fitted. This method yields a factor at all posi- 
tion angles, which is nominally equal to 1 for the subset of data 
at PA 70°, and 1.16+0.01 and 1.23±0.03 for the data along PA 
170° and 130°. Assuming that the circumstellar structure around 
MWC 297 is a geometrically flat disk - supported by the low 
L IR / L* ratio of 5% - these numbers can be translated into an 
inclination and position angle of the major disk axis on the sky 
of 40+10° and 120+50° respectively. The extent of the com- 
ponents along the disk major axis is 1.3+0. 1 times larger than 
along PA=70°. The reduced chi-square of the rescaled model 
visibilities, compensated for the additional model parameter, is 
3.8 (bottom part Table [5]). Although a large inclination is ex- 
cludable, we caution against overinterpretation of the inclina- 
tion and disk position angle estimates. The two baselines with 
different PAs are also the longest baselines, and hence sample 
a slightly different spatial-frequency domain than the subset. A 
perfectly spherically symmetric model may therefore fit the data 
just as well (as was e.g. shown with the MIDI fits). To straighten 
out the issue and determine the exact inclination, more visibility 
measurements at comparable spatial frequencies, but different 
baseline angles, are required. 

We consider our fit to the subset of data as our most reli- 
able fit, since no issues related to the position angle are present. 
The model fit to the visibility data is shown in Fig. [4] In Fig. [5] 
a radial cut of the intensity distribution of the model at differ- 
ent wavelengths is shown. Note that, although the FWHM of the 
largest GD component is 40 mas (10 AU), the FWHM of the to- 
tal image is less than 6 mas (1.5 AU) at all wavelengths. The ex- 
tent of our approximated image of MWC 297 is in perfect agree- 
ment with the VISIR image at 11.85 /mi. We have convolved 
the model image with a Gaussian curve of FWHM 0.33±0.01", 
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Fig. 4. Top panels: Best model fit to the visibility data. The AMBER data are corrected for the stellar flux contribution. The two 
longest MIDI baselines have different baseline position angles. The model fit to the visibilities related to the latter is the best-fit 
model which includes the effect of the disk's inclination, as described in the text. Bottom panels: Residual difference between 
observed and modeled visibilities, expressed in terms of the error cr on the individual measurements. [This figure is available in 
color in the electronic version.] 



to mimic the point spread function estimated from the two cal- 
ibrator measurements. The resulting image has a FWHM of 
0.39+0.01", in perfect agreement with the measured 0.41". 

We have compared the H-band AMBER a nd model visi- 
bilities to the broad-band visibilities obtained bv lMonnieret alj 
(2006) with IOTA. In Fig.|6l their H-band visibilities are shown 
together with the AMBER visibilities at 1.65 //m and the broad- 
band visibilities of our best-fit model. The IOTA measurements 
are in relatively good agreement with our data given the large 
error bars, although they appear to be systematically higher. It 
is not clear to us why this discrepancy is present. The IOTA H- 
band data were obtained at different baseline position angles. 
The smooth decline of the visibilities with spatial frequency, in- 
dependent of position angle, shows that the circumstellar emis- 
sion region does not display strong asymmetry. In case of a flat 
disk, this corresponds to a moderate inclination of at most 40°, 
consistent with our upper limit. 

3.3. Consistency with the SED 

The brightness temperatures derived for the three components 
in our best-fit model, are consistent with the near- to mid-IR 



SED. The model emission accounts for the total near- to mid- 
IR flux. In Fig. [7] the contribution of the components to the 
total flux is shown. The spectrum of the largest component is 
a single-temperature Gaussian disk consisting of olivines with 
a column density of 0.02 g cirT 2 . At sub-mm/mm wavelengths, 
only the coolest component contributes significantly. The pho- 
tometric points at these wavelengths are clearly higher than the 
predicted model flux. This points to the presence of an additional 
cold dust component, which does not contribute flu x at near- to 
mid-IR wavelengths. Following Manoi et al. (2007), 200 mJy of 
the 1.3 mm flux can be attributed to dust emission, while the 
rest is free-free emission from an ionized wind. After subtrac- 
tion of the contribution of the three warm near-/mid-IR compo- 
nents, 155 mJy remains. The combination of the 1.3 mm flux 
and the fact that the cold dust does not contribute to the mid- 
IR flux, places an upper limit of 200K on the brightness tem- 
perature of this material, and a lower limit of 160 mas to its 
extent. The limits set by spatially resolved mm observations 
(< 230 mas, lManningsiri994t iManoi et aJ1l2007l) . constrain the 
brightness temperature to more than 100K. The brightness tem- 
perature of the cold dust is therefore restricted to the interval 
between 100 and 200K. 
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Spitzer IRAC images of MWC 297 show very complex cloud 
structure on arcminute scales (Fig. [8]). Its extent and structure 
coincides with that of the 50-1 00 /mi brightness distribution 
measured by di Fr ancesco et al.l dl994l) . This cold material is 
likely responsible for the additional far-IR (>25 /mi) flux ex- 
cess, evident in e.g. the IRAS 25, 60 and 100 micron points. Also 
the ISO-SWS and ISO-LWS (Long Wavelength Spectrograph) 
spectra, which are obtained with increasin g aperture size to- 
wards longer wavelengths dBenedettini et alj|2001l) . suffer from 
this large-scale emission. It is possible that the mm observations 
are partially affected by this emission as well, although most of 
it would be resolved out. At near- and mid-IR wavelengths, the 
flux contribution of the clouds is negligible: at 8 micron, the to- 
tal flux adds up to 130 Jy, spread out over roughly 10 square 
arcminutes. Assuming a homogeneous emission region, the flux 



inside the ISO 14" x 20" aperture would be of the order of 1 Jy 

(<1%). 

4. Discussion 

4.1. The extinction towards the MWC 297 system 

The flux of the central star in th e MWC 29 7 syste m is at- 
tenuated with 8 magnitudes in V. iDrew et all (1997) suggest 
that the real extinction could be even stronger, if the observed 
stellar flux is scattered, rather than direct light. However, this 
seems unlikely; The application of an interstellar extin ction law 
(ISavaee & Mamislfl979l: ISteenman & Thdfl98l fl99ll with A v 
= 7.75 and Ry = 3.1), makes that the dereddened photometry 
agrees with the expected fluxes of a B1.5V star at 250 pc from 
UV to the optical. In the case of scattering, the observed pho- 
tometry would have been bluer. Application of an interstellar 
extinction law would in this case be inappropriate and would 
yield results inconsistent with a B1.5V model. 

A priori, it is unclear whether the attenuation of the central 
star is due to circumstellar matter, physically linked to the cen- 
tral star, or interstellar dust. However, in the case of MWC 297, 
the extinction is most likely interstellar. An extended CO cloud 
is detected in the direction of the target (Manoj et al. 2007). The 
extent of this cloud is at least 1' x 1', much larger than the ex- 
tent of the near- and mid-IR emission region of MWC 297. The 
ISO spectrum of MWC 297 shows H 2 and C0 2 ice absorp- 
tion bands, indicative of a cold environment. It is improbable 
that these ices are abundant in a small (R < 80 AU) and hence 
warm environment close to a B1.5V star. Furthermore, attenua- 
tion of the central star can be descibed with an interstellar ex- 
tinction law, where small (0.1 /mi) grains are the main source 
of opacity. The presence of small dust grains, CO2 and H2O 
ice, and the large-scale CO emission are naturally explained by 
an interstellar molecular cloud in between the observer and the 
MWC 297 system. Moreover, the proximity (0. 1 pc at 250 pc) 
of the large-scale cloud structure present in the Spitzer images 
-which is als o evident from the exte nded 50-100 //m emis- 
sion found bv ldi Francesco et al] ( 1 19941) - strongly suggests that 
the latter is remnant material from the natal cloud f rom which 
MWC 297 and the nearby lower-mass young stars dVink et al.l 
2005; Damiani et al. 2006) were formed. It is probable that the 
interstellar cloud in the line-of-sight is part of this large structure 
surrounding the MWC 297 stellar group. 

The main consequence of this conclusion for the present pa- 
per, is the fact that the entire MWC 297 system (i.e. the central 
star + the circumstellar matter responsible for the near- to mid- 
IR flux) is screened by this cloud. This is not important for the 
interpretation of the interferometric data, because the latter are 
only sensitive to the normalized brightness distribution on the 
sky. Relative quantities (e.g. flux ratios between two components 
in the interferometric models) can however be transformed into 
absolute quantities (e.g. brightness temperatures), only if the un- 
obscured spectrum of the system is known. We have dereddened 
the AMBER, MIDI and ISO spectra of MWC 297 for that pur- 
pose, with the same extinction law extending from UV to mid- 
IR. 



4.2. The near- and mid-IR geometry of MWC 297 

Our investigation of the near- and mid-IR brightness structure of 
the MWC 297 system yields the following results: 

- MWC 297 is an extremely compact source. The bulk of the 
near- to mid-IR emission emanates from a region remark- 
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ably close to the central B1.5V star. The extent of the sys- 
tem in the H-, K- and N-bands (image FWHM ~ 1.5 AU) 
is well within the maximal outer radius of t he system deter- 
mined from mm ob servations (< 60-80 AU. lManningsll994t 
iManoi et alj2007l) . 

The closure and differential phases obtained with respec- 
tively AMBER and MIDI are equal to zero within the er- 
rors. There is no evidence of strong deviations from centro- 
symmetry along a position angle of 70° E of N at the spatial 
resolution of our observations, neither in the continuum nor 
the silicate feature. 

There is no significant IR-emission-free gap between the cir- 
cumstellar matter and the stellar surface at the resolution of 
our observations. 

Tentative evidence is found for a moderate elongation of the 
brightness distribution, or a moderate inclination of the cir- 
cumstellar disk (z < 40°). A high inclination is excludable. 
Parametrized accretion disk models cannot explain the data. 
This either suggests that the circumstellar matter around 
MWC 297 does not reside in an accretion disk, or that sim- 
ple models are not appropriate to explain the characteristics 
at the small angular scales observed. 



- A spherically symmetric dusty outflow or halo cannot ac- 
count for the interferometric and photometric observations 
either. 

We have used geometrical models to approximate the cir- 
cumstellar brightness distribution of MWC 297, with as few 
free parameters and analytic components as necessary. Three cir- 
cumstellar components suffice. However, it is unclear what the 
physical meaning of these components is. Although it might be 
that they are independent physical structures, it is more likely 
that the three Gaussian disks sample the temperature gradient 
of a smooth density distribution. Our analysis has shown that 
a power-law approximation of the temperature is inappropri- 
ate for the circumstellar matter around MWC 297. Furthermore, 
the strong radial drop-off of the temperature seems to suggest 
that the circumstellar matter is optically thick in the radial di- 
rection. The halo modeling has shown that an optically thin or 
moderately thick circumstellar medium produces a much shal- 
lower temperature gradient: at the location of the largest geo- 
metrical component, the temperature would be much higher than 
observed. The most likely geometry of the circumstellar matter 
is therefore a flattened disk, which is viewed under a moder- 
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ate inclination. However, this disk is significantly different from 
a typical lower-mass Herbig Ae or T Tauri accretion disk: it is 
compact and geometrically rather flat (as supported by the low 
IR excess luminosity). 

4.3. On the nature of the circumstellar matter 

The circumstellar structure in the MWC 297 system is remark- 
able in many respects. We highlight some literature results and 
discuss these in the framework of the present paper. 

- The central star rotates with a pro jected velocity vsin/ = 
350 ± 50 km s -1 dDrewet alii 19971) . Critical rotation occurs 
when v rot attains yflGMJOR^) = 450 + 70 km s" 1 , or 
equivalently when i < 55 ± 15°. The interferometric mea- 
surements suggest a moderate inclination i < 40 + 10°, and 
hence that the star is indeed rotating at or close to critical 
velocity. 

- Extremely strong hydrogen emission lines from optical to 
mid-IR wavelengths are observed. The spectral line profiles 
are rather narrow (Ho- FWHM ~ 200 km s _1 ) compared to 
e.g. the half-maximum widths observed in other Herbig Be 
stars. The ex tinction of the hyd rogen lines is consistent with 
an Ay of 8 dDrew et al.lll997l) . which suggests that the hy- 
drogen emission region is attenuated by the same material as 
the central star. 

- The forbidden [Oi] emission lines at 6300 and 6363 A 
are extremely str ong, but al s o extr emely narrow (FWHM 
~ 22 km s" 1 , lAcke et all 120051) . The lines are ther- 
mally excited, as proven by the detection of the [Oi] 
5577 A line (lAndrillat & Jaschekl [l998h : The line ratio of 
the 6300 and 5577A lines is 5+1, combining the equiv- 
alent width measurements of lAndrillat & Jaschekl dl998l) 
and Acke et al. (2005). Typical Herbig Be so urces (e.g. 
Z CM a,CD-42°11721, PV Cep, V645 Cyg; lAcke et all 
2005) display strong [Oi] lines with extended blue wings, 
up to several hundreds of km s -1 . The latter is a strong indi- 
cation for the presence of a bipolar outflow. If the [O i] emis- 
sion in MWC 297 also emanates from a n outflow, the outf low 
axis should lie in the plane of the sky (Ac ke et al.ll2005l) . to 
match the limited velocity dispersion of the observed spec- 
tral lines. Given the small disk inclination, however, it is 
most probable that the [O i] emission emanates from the disk 
itself. 

- iManoi et al.l (120071) find a shallow sub-mm slope of the SED, 
indicative for the presence of large grains, or an emission 
region which is optically thick at these wavelengths. The first 
would indicate grain growth, the former is consistent with a 
very compact and dense circumstellar disk. 

- In the 3-micron range, the spectra of MWC 297 display dia- 
mond emission dTerada et al.ll200lt lAcke & van den Anckerl 
2006). Diamonds form under conditions of severe tempera- 
ture and pressure, both of which appear to be present in the 
circumstellar environment of MWC 297. 

- H2O maser emission has been found by lHan et al.l y998) 
at a Vlsr of -75.6 km s -1 . The emission is blueshifted by 
~ 85 km s compared to the radial velocity of the central 
star (vlsr = 10 ± 2 km s -1 , based on t he pea k position of 
optical emission lines). Palla & Prusti (1993) suggest that 
in Herbig Be stars the origin of this emission lies in a col- 
limated (bipolar) outflow. The detection of water emission 
seems to suggest that also MWC 297 drives a molecular out- 
flow, despite the absence of a detected CO outflow, even at 
high blueshifts (see Fig. |3). Given the strong blueshift, the 



outflow axis of the H2O must be approximately aligned with 
the line of sight. The receeding (redshifted) emission would 
then be obscured by the circumstellar disk. This is consistent 
with a small disk inclination. 

MWC 297 is a high-mass young star, accompanied by lower- 
mass newborn stars and remnants of its natal cloud. The cen- 
tral B1.5V star is a critical rotator. It is surrounded by a dense 
gaseous and dust-rich disk, which reaches in close to the stel- 
lar surface and produces the prominent gaseous emission lines 
in the optical-to-IR spectrum. The star+disk are viewed under a 
moderate inclination angle. 



5. MWC 297: a critical rotator and Be star 
progenitor 

Stellar rotation has, in recent years, been realised to be a 
fundamental property that can deter mine the evolution and 
ultimate fate of high- mass stars dMevnet & Maeden [2000; 
Maeder & Meynet 2000). Rapid rotation impacts on stellar life 
times and alters the chemical evolution of the star. Knowing 
stellar rotation as a function of mass and age is thus essential 
to understanding high-mass star evolution. In our galaxy, some 
twenty per cent of all high-mass stars, the classical Be stars, 
show rapid rotation and associated mass loss. Our identification 
of MWC 297 as a critical rotator proves that high-mass stars can 
achieve such high rotation rates at birth. In this section, we dis- 
cuss the consequences of this discovery. Furthermore, we argue 
that MWC 297 is an evolutionary progenitor of a classical Be 
star. 

Observations have provided increasing evidence that also 
high-mass stars form due to mass gain via an accretion disk 
dVink et alj|2002l) . as do solar-mass stars. For the latter, the rota- 
tional velocity at birth is thought to be determine d by the interac- 
tion of the star and its accretion disk (S hu et alj|1994l) . Magnetic 
fields lock the stellar surface rotation rate on the orbital rate close 
to the inner edge of the disk. A high accretion rate, as well as 
a low magnetic field strength, shifts the distance at which stel- 
lar and disk rotation are locked to smaller radius, implicating a 
higher stellar rotation rate. Therefore, theory predicts that stars 
with a small inner disk radius have a high surface rotation rate at 
birth. 

The evolution of the equatorial rotation speed of high-mass 
stars is governed by the combined effects of surface mass loss 
and the internal redistribution of angular momentum. The first 
removes angular momentum and is believed to be important for 
very massive (M>25 M Q ) stars; its influence for lower-mass stars 

1 i ; ^\\ 1 

is unclear (Ekstrom et al. 2008) but can be significant. Angular 
momentum redistribution inside the star, however, can result in 
an increase of the ratio between equatorial rotation speed and the 
critical velocity. Theory predicts that the critical velocity may 
be attained for stars of 10-20 M at a considerable age of 10 
million years. On the other hand, very little is known about the 
rotational velocities of high-mass stars at birth, and theory is cur- 
rently unable to describe the evolution of stars born at the crit- 
ical speed. It is in this framework that the critical rotation of 
MWC 297 attains its full significance. 

Two observations imply that a dense and opaque gaseous in- 
nermost disk, very close to the stellar surface, shields the outer 
disk in the MWC 297 system. (1) The occurrence of dust so close 
to the star is remarkable, as direct stellar irradiation would heat it 
there far above its sublimation temperature. In lower-mass young 
stars, stellar irradiation makes the innermo st disk dust-free an d 
opens up a gap between star and dust disk (Monnier et al. 2005). 
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Fig. 8. Traveling to the center of MWC 297. Left: Spitzer IRAC image at 8 /im. Clearly visible is a complex, large-scale cloud 
structure. A numerical circular mask covers the saturated part of the image. Middle: Zoom-in at the center of the Spitzer image: 
VISIR image at 11.85 /xm. The image is seeing-limited, but its FWHM is significantly larger than that of the calibrators. Right: 
Zoom-in at the center of the VISIR image: our best-fit interferometric image at 12 yum. The grey scale in the three images is linear. 



No such gap is observed in the disk of MWC 297. (2) Given 
the strength of the optical emission lines but their small radial- 
velocity dispersion, these lines must mainly originate in the disk, 
which is at low inclination. A large amount of neutral oxygen 
and hydrogen atoms are thus present close to the star, which 
confirms that the ionising stellar photons are absorbed efficiently 
close to the star. 

The fact that this remarkably compact disk encompasses a 
critically rotating star suggests a causal link. This link could be 
the formation scenario discussed above. In that case, MWC 297 
would either have a weak magnetic field, have experienced high- 
rate mass accretion, or a combination of both. It would mean that 
a scenario initially developed for low-mass stars, also applies for 
their high-mass counterparts. Accretion may still be ongoing in 
this case. The current observations, however, do not allow for a 
definite judgement on whether the matter in the innermost disk 
is moving inward or outward. It cannot be excluded that accre- 
tion has ceased, and that the critical rotation of the star currently 
induces mass loss to the disk through an equatorial outflow. Both 
hypotheses are consistent with the absence of an inner gap im- 
plied by the near- and mid-infrared observations, and the pres- 
ence of high-density gas very close to the stellar surface. 

The discovery of a young star rotating at critical velocity 
proves that stars can achieve such high rotation rates at birth. 
This conclusion naturally raises intriguing questions on the ini- 
tial distribution of rotational rates in newly formed high-mass 
stars. At this point, it is difficult to estimate the frequency at 
which critically spinning stars are produced. A hint that this 
fraction cannot be too small may be found in the proximity of 
MWC 297: it is the closest high-mass young star known to date, 
and as such the prime target for investigation. That this partic- 
ular star rotates at critical velocity suggests that we were either 
very fortunate, or that fast rotation at birth is more common than 
generally assumed. 

It is tempting to speculate about the future evolution of 
MWC 297. The star will most likely evolve into a classical Be 
star. The latter are high-mass main-sequence stars which ex- 
hibit line emission produced in a gaseous equatorial disk: the 
line profiles of the Balmer lines indicate Keplerian rotation and 
rather small outflow velocities. Since classical Be stars share the 
common characteristic of faster-than-average stellar rotation, it 
is likely that strong equatorial centrifugal forces contribute to the 



origin of the disks. However, observed rotation velocities do not 
reach the critical velocity, hence the need to invoke additional 
effects, su ch as puls ations or magnetic streams, to explain disk 
formation dPorter & Rivin ius 2003). In a first phase, MWC 297's 
disk will evolve from dust-rich to purely gaseous: The high den- 
sity in the innermost disk will decrease, either due to a decrease 
of the accretion rate, or a decrease of the stellar equatorial mass 
loss. The disk will eventually become transparent to stellar pho- 
tons, which will sublimate the dust grains. The neutral oxygen 
and hydrogen will become ionised, implying the disappearance 
of the forbidden lines and significantly reducing the strength of 
hydrogen emission. MWC 297 will then have all the features of 
a classical Be star. 

Our study demonstrates the power of interferometry to de- 
termine the initial rotation rates of massive stars. When such 
studies can be extended to a statistically significant sample, it 
will become possible to assess whether MCW 297 is an excep- 
tional object or whether critical rotation is the general rule for 
newborn massive stars. In the latter case, the much slower ro- 
tation of most main-sequence stars would imply that the pro- 
cesses which carry away angular momentum are very efficient 
in our galaxy. In the former case, MWC 297 would qualify as 
a candidate for the hypothetical objects in which extreme rota- 
tional mixing induces an almost homogeneous internal composi- 
tion throughout the main-sequence pha se, hence c ausing the star 
to bypass the red-supergiant phase dMaederill9~87l) . MCW 297 
would then be a possib le progenitor of a gamma-ray burster 
dWooslev & Hegeri2 006) close to our sun. However, gamma-ray 
bursters appear to occur preferentially in low-metallicity envi- 
ronments (Wolf & Podsiadlowski 2007), where indeed the mass 
loss and associated angular-mome ntum losses are expe cted to be 
lower than in our galaxy (Mac Fadven & Wooslevl[l999l) . 
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